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Abstract 

The fusion cross section of 3 He( 3 He,2p) 4 He at a center of mass energy of 30 to 50 keV has been measured 
by using helium-3 doubly ionized beam at a low-energy high current accelerator facility, OCEAN. Free 
from molecular interference in the beam, the measurement determines the astrophysical S-factor with better 
statistical and systematical errors than previous data. By using singly and doubly charged helium-3 ions the 
facility envisages to provide the data from high energy to Gamow energy regions. 
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I. INTRODUCTION 



Of the reactions involved the solar combustion, namely d + p — > 3 He + 7, 3 He + 3 He — > 2p 
+ a, and 3 He + a — > 7 Be + 7, we have focused on the cross section measurement of the 3 He + 
3 He reaction at an effective energy of E cm =30-50 keV. Among many experimental works, the only 
experiment which has been conducted in or lower than the Gamow energy is the 3 He + 3 He — > 
2p + a reaction at the LUNA in Gran Sasso laboratory [1]. The latter was performed near the 
Gamow energy without measurement at E cm =25-45 keV due to accelerator limitations that lower 
voltage could not be applied with a high voltage static accelerator while higher voltage more than 
50 kV could be impossible at the LUNA. For a wider energy range E cm =17.9 to 342.5 keV, the 
experiments done by Krauss et al slightly extended the data over this energy gap |2D. For the 
nuclear astrophysics discussion, in particular, standard solar model or nucleo-systhesis continuous 
data down to low energy is crucial to deduce the astrophysical S-factor, S17. Therefore, successive 
and precise data from around 50 keV to 20 keV center of mass energy are needed. For this purpose, 
we have constructed a low energy and compact accelerator facility which provides doubly charged 
3 He ions for the measurement in the region between 50 keV to 25 keV and also singly charged 
3 He ions in the region less than 25 keV center of mass energy. This is the first report of a series 
experiments done at our low-energy high current accelerator facility, OCEAN. We obtained precise 
results from 45.3 to 31.2 keV in center of mass energy for the 3 He+ 3 He^2p+a reaction with a 
doubly charged incident 3 He beam. It has the substantial advantage of eliminating the molecular 
interference in the proton spectra by the 3 He + d reaction caused by the incident HD + beam. 

II. EXPERIMENTAL APPARATUS 

The experimental apparatus, OCEAN (Osaka University Cosmological Experimental Appara- 
tus for Nuclear Physics) consists of: (1) a powerful ion source that provides an intense current of 
more than 1 mA for 3 He 1+ ions at incident energies of 30-50 keV, or 200 fiA for 3 He 2+ at incident 
energies of 90.6 to 62.4 keV. (2) low energy beam transport with good transmission, (3) a win- 
dowless gas target and recirculation / purification system, (4) a reliable calorimeter, (5) detectors 
for reaction identification, and (6) an electronics and data acquisition system based on CAMAC. 
The layout of the OCEAN is shown in Fig. [TJ 
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A. Ion source and extraction electrodes 



An intense ion source that can produce 3 He 2+ ions is essential for the present study. The 
NANOGUN™, which was obtained from PANTECH, confines high-temperature electrons pro- 
duced by the electron cyclotron resonance (ECR) and is assembled into an ECR ion source with 
10 GHz, 200 watt RF generator (model VZX-6383G5, CPI). The original NANOGUN ECR ion 
source made by PANTECH can provide 40 fiA for 40 Ar 8+ at 20 kV extraction voltage with an 
RF power of 60 W. From these data we could easily foresee the possibilities of obtaining 3 He 1+ 
or 3 He 2+ ion beams of more than 100 fx A, although an applied potential to the extraction is not 
enough to use this ion source for astrophysical applications in a wider enery range. For this pur- 
pose, we redesigned the ion source extraction system to meet the ion optical condition for the 
present window-less gas target. Leroy et al. reported an improvement for the analyzed current and 
beam emittance of 3 He 1+ beam by the a system called Multielectrodes Extraction for CAPRICE- 
type ECR ion source. It supplied with an analysed current of 2.3 mA of 3 He 1+ giving a transmis- 
sion of more than 75 % in the beam line [3]. did you also do this). Based on several experimental 
studies and computer simulations we designed and applied a two-electrode extraction system for 
the NANOGUN ECR ion source [4]. This improves the beam emittance under the influence of a 
strong space charge force, and secondly, it moderates the electric field gradient ascribed to high 
operational voltage. The optimization process for applying a suitable voltage to the intermediate 
electrode progressed considerably during experiments measured at each center of mass energy. 
For the present experiment between 45 keV and 31.2 keV center of mass energy, around 100 /iA 
of double charged 3 He beam could pass through the three stages collimators to the gas target. A 
specially designed ceramic insulator (SUMIKIN CERAMIC Co-Ltd) with two folds on the surface 
could support up to 50 kV potential without any surface discharge at a distance of 170 mm from 
the surface. We fabricated conic type and straight type extraction electrodes. It consists of a beam 
forming electrode which has a long straight part similar to the original NANOGUN's design and 
a conic structure for the beam forming electrode nearly parallel to the end surface of the plasma 
chamber. The performance and design have been described in previous papers \^\. 
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B. Low energy beam transport 

The low energy beam transport system between the ECR ion source and the gas target achieved 
a high beam intensity and other desirable beam qualities, thereby allowing for precision measure- 
ments of the present 3 He+ 3 He— >2p+a experiment. Generally, it is known that there a strong space 
charge effect in the beam transport at ion currents of more than 1 mA. It is essential that this effect 
is accounted for when calculating the beam optics. We used a GIOS code developed by Wollnick 
et al. for incorporating this effect yfl. We adopted the D (dipole, 90-degree deflection angle) + Q 
+ Q transport scheme for our system since it is easier to operate fewer elements. In this calcula- 
tion we assumed a beam source of 100 n mm mrad and 5mm in diameter. Despite the variance 
in total potential of the beam, a nearly invariant beam form could be realized at the source exit 
using the extraction system stated before. To maintain the minimum slit aperture, we calculated 
the dimension of the beam at the target position by varying the various parameters of elements and 
drift lengths so as to achieve smaller dx and dy. Very attractive results such as a constant dx and 
dy and a nearly parallel beam form are shown in Fig. El The beam transmission efficiency from 
the ion source through the target is about 30 % (Table[Q). 

C. Window-less gas target 

The window-less gas target for the study of the 3 He+ 3 He reaction consists of a differential 
pumping and gas circulation/purification system as shown in Fig. |3] 

In order to maintain a pressure of 1 mbar in the chamber of the 3 He gas target without window, 
the pumping system should be composed of several stages between the target chamber and beam 
transport system. Thus we prepared a helical groove vacuum pump (model TS-440, OSAKA 
VACUUM Co. Ltd TS-440) as the main pump for evacuating the gas flow at the viscose region as 
well as at the higher vacuum region. The sizes of collimators at every stage were estimated by 
calculation in order to maintain a pressure of about 0.3 Torr in the target chamber. During the 
experiment, the pressure in the target chamber was maintained at 3.1(±0.1) Torr for about four 
days. 

The recirculation system consists of a helium tight pump, an oil free diaphragm membrane 
compressor, reserver vessels, compound gauges, ultra fine regulated valves and a quadruple mass 
spectrometer as shown in Fig. Since all evacuation pumps connected to the helical pump, a gas 
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dosing system has been employed in order to maintain a constant target pressure to evacuate the 
input gas amount at the ion source. Because during experiments, gas injection into the ECR ion 
source had the effect of gradually increasing the target pressure. Thus, a constant target pressure 
was maintained by controlling the signal from the capacitance manometer (Barocel-655) located 
at the target chamber with a gas dosing apparatus (model EVR 116 and RVC 300 controller). 

The new purification system developed for the present experiment is quite different from the 
usual method. It exploits a cryo-pump (model U-140W, DAIKIN) as a purifier without a special 
absorbent and liquid nitrogen, although high heat input could be expected at the high operat- 
ing pressure. This has been overcome by adding another oil free turbo-molecular pump (model 
A30FC, ANELVA) between the cryo-pump and the target chamber. As pointed out by A. Krauss 
et al. [2], the deuterium contamination both in the target and the beam, resulting from the water 
vapor, is a crucial problem for obtaining low energy data, since the d + 3 He reaction cross sec- 
tion is six orders of magnitude higher than that of the 3 He + 3 He reaction. In the case of 3 He 2+ 
beam (e/m=2/3), we can avoid molecular interference of the HD + beam (e/m=l/3), with analyzer 
magnet (Fig. [TJ contaminating the incident 3 He 1+ (e/m=l/3) beam, whereas this could not be 
avoided if we had employed a single charged 3 He beam (e/m=l/3), we have measured the deu- 
terium contamination in commercial 3 He gas by detecting HD + separately via Accelerator Mass 
Spectrometry (AMS). The experiment was carried out using the RCNP K=140 AVF cyclotron. 
The cyclotron accelerator and the NEOMAFIOS ECR (NEOMAFIOS) ion source were operated 
only for the experiment on the beam injection line at an anex leading to the post accelerator (RCNP 
Ring cyclotron). The present result is HD + / 3 He = (3.82±0.69)x 10~ 5 . It is noticed that even if an 
ECR ion source uses electrons with higher energy than that of a duo-plasmatron ion source there 
is a considerable amount of HD + production from the surface of the ion source and from the bottle 
of 3 He gas. 

Deuteron contamination in the target was also estimated during operation of the recirculation 
and purification system. Then the pressure at the target chamber was 1 .2 x 10 -2 Torr, and 7.6 x 10~ 7 
torr at the helical pump. The H 2 component in the residual gas was measured by means of 
a quadrupole mass spectrometer, to be about 20 %. Assuming that the amount of H 2 at the 
target gas is the same as the residual gas, and the deuteron abundance is the same as the natural 
abundance (0.014 %), we can deduce that the deuteron contamination (D 2 0) is in the order of 
ppm. This is satisfactory for the present measurement as will discussed later. In addition, the 
amount of deuterium contamination in the target gas was also evaluated by means of the 3 He + 
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3 He experiment, by detecting the 14.7 MeV proton, and was found to be about 0.1 ppm. 
D. Calorimeter 

A calorimetric device has been developed to measure the projectile flux in the present experi- 
ment, since the conventional method of charge integrating to determine the number of incident par- 
ticles is not applicable due to the neutralization of the incident charged particles with the gaseous 
target. 

Many types of calorimeters have been developed There are two types of calorimeters. One 
type measures the temperature difference between two parts thermally isolated with large heat 
resistance material. The other type measures the power needed to stabilize the temperature of 
water used to cool the calorimeter cup which is bombarded with energetic particles. 

We developed the calorimeter by using a heat flux sensor with an accuracy of better than 5 % 
in the dynamical range of 1-30 watts [8]. 

The structure of calorimeter is shown in Fig. |U It consists of a solid copper heat sink (100 
mm length, 38 mm diameter) with water channels, and a Faraday cup (140 mm length, 38 mm 
diameter, 1.5mm wall thickness) in front of the heat sink. The Faraday cup is supported by a 
Pyrex glass insulator as well as a metal flange with a Teflon gasket. This organic gasket served 
as a vacuum seal as well as an electric insulator. Therefore, this calorimeter can also measure the 
number of incident particles, when the target chamber is in vacuum, by an electrical method. At 
the entrance of the cup, a secondary electron suppressor is installed. Around 100 volt was applied 
to the suppressor for the current measurement. 

After passing through the window-less gas target, the incident 3 He 1+ and 3 He 2+ beam is neu- 
tralized and captured in a Faraday cup and calorimeter. In order to be applicable to a wide current 
range, a heat flux sensor (HFS) (OMEGA HFS-3) was used to measure the heat transfer from 
the hot part to the cold part. The HFS is a thermister bolometer which can measure the heat 
flux to or from a surface with minimum disturbance of the exiting heat flow. In this method, a 
self-generating thermopile is arranged around a thin thermal barrier to produce a voltage that is 
a function of the thermal energy passing through the sensor. The response of the sensor to the 
thermal energy input is 1.10 to 1.11 mVAV/m 2 . For precise measurements of the heat transfer, we 
made use of the following procedures: 1) In vacuum, the current of the ion beam can be measured 
by standard charge integration where the calorimeter cup serves as a Faraday cup. 2) To reduce 
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the conduction and convection losses, the heat capacity of the calorimeter should be reduced by as 
much as possible in order to get a better time response. With the present heat flux sensor, tempera- 
ture differences of less than 0.001 °C can be detected easily. Therefore, the temperature of the heat 
sink of the calorimeter should be stabilized to better than 0. 1 °C with heat exchange. As shown 
in Fig. HJ between the copper base and the thin plate of the calorimeter, two heat flux sensors are 
sandwiched with thermister temperature sensors. These are originally insulated electrically. As 
shown in Fig. |4j termister temperature sensors are also located to measure the temperature of the 
ambient or Faraday cup base and the cooled heat sink. These are installed in a stainless steel pipe 
(40 mm length, 10.5 cm diameter), which can be evacuated by a small turbo-molecular pump. 

The calorimeter was tested by using a 3 He 2+ beam of energy of 40 kV. The relation between 
beam current(I) and heat flux(H) can be written as, 

I ■8t = k 1 - H -6t + C -ST, (1) 

where T is the temperature of the calorimeter and C is the heat capacitance. The term C ■ 5T 
shows that the temperature of the calorimeter depends on the incident beam current. Thus, if the 
intensity I is changed, the converted heat is used to heat the calorimeter and is also transferred to 
the cold base. 

Since the transferred heat may approximated to be proportinal to the temperature difference 
between the the front and cold base, the second term in equation [l] C ■ 5T can be rewritten as 
k 2 ■ 5H. Thus, equation (1) can be written as, 

I = h-H + k 2 S -^. (2) 

In order to determine the parameters k\ and k 2 , an experiment was carried out with a 3 He 2+ 
beam of 40 keV ( E cm ). The beam current was calculated from the HFS output by comparison 
with the Faraday cup measurements. The HFS output was measured by a KEYTHLEY 2000 
multimeter. The parameters k\ and k 2 were determined as follows: 

1) Parameter ki : If the system is stable, that is, incident beam / and temperature T are stable, the 
second term of Eq. (2) can be ignored. In this condition, the parameter k\ can be determined by a 
least square fit as shown in Fig. |3a). 

2) Parameter k 2 : This papameter can be determined once parameter k\ is known, as seen in Eq. 
(2). The term 5H/5t was measured for averaged time scales of 3, 7, 15, and 30 s. It was found 
that the scale of 30 s was appropriate. Fig. |5tb) shows 8H/8t as a function of / — k x ■ H. 
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3) Comparison with beam current: Fig. |6] shows the beam current as a function of time measured 
using a Faraday cup, k\ ■ H, and k\ • H + &2^§- with different beam conditions, i.e., stable, slowly 
increasing, and decreasing beam current. 

The calculated currents measured with the HFS reproduce the measurements with the Faraday 
cup as shown in Figs. EJa) to|(Jc). On the other hand, if the beam current is suddenly changed as 
shown in Fig. Eld), the calculated currents from the HFS output overestimate the value measured 
by the HFS. Further improvements of this system are necessary. Fortunately, this case should not 
be a fatal problem for an astrophysical experiment with long term measurement times of typically 
one day or one month, since this occurs rarely, at most one or two times in a day. Thus, in the 
experiment, the error caused by this overestimate can be neglected. 

1. Reproducibility of beam current for different energies 

The reproducibility of the beam current determined from the HFS output was verified for sev- 
eral beam energies. Measurements were carried out on 3 He 2+ beam at incident kinetic energies 
of 35, 30, 25, and 20 keV. The beam currents were calibrated using the parameters determined 
with the 40 keV beam. Fig. |7] shows the accuracy of the calculated beam current in the form of 
(Ihfs — Ifc)/Ifc, where Ihfs and Ifc are the beam currents measured by the HFS and the 
Faraday cup (FC), respectively. The accuracy was better than 2 % for the measured energies. 

2. Estimation for transferred heat in HFS 

The heat transferred through HFS can be calculated from the calibration parameter of HFS 
(OMEGA HFS 3: 1.10 ~ 1.11 /iV/W/m 2 at 70 °F). Also the kinetic energy of incident particles 
can be calculated from the charge integration. Table Hll shows the results of these calculations. It 
was found that in the equilibrium state most of the heat is transferred through the HFS. 

3. Heat exchange with the surroundings 

Because of the different vacuum conditions during the calibration (~10~ 6 Torr) and during 
the experiment (0.1 Torr), the effects of convection heat losses by the flowing gas in the target 
chamber have to be estimated. The heat transmission by convection was measured by comparing 
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the transferred heat through the HFS for the two conditions of vacuum (10 Torr and 0.1 Torr) in 
the target chamber. The temperature of cool base was 30 °C. The difference of HFS output was 
0.04 W corresponding to about 1.7 % for a 35 keV 100 fiA beam. 

E. dE-E counter telescope 

In order to ensure a large detection efficiency and a clear discrimination of real events, we 
exploit four dE-E counter telescopes by using semiconductor detectors for the measurement of the 
3 He( 3 He, 2p) 4 He reaction. These detectors are installed into the target chamber filled with 3 He 
gas and are capable of identifying the 3 He+ 3 He reaction (Q = 12.86 MeV) as shown in Fig. [8] 

The reaction generates two protons which have kinetic energies of to 10.7 MeV, and an alpha 
particle which has kinetic energies of to 4.3 MeV. The dE and E detectors in each telescope 
have an active area 2500 mm 2 , the dE detector has a thickness of 140 /im and the E detector has a 
thickness of 1500 /mi (MICRON Ltd. ). 

To stop the generated alpha, photons and elastically scattered 3 He from the beam, an aluminized 
Mylar films with thickness of 25 /im are located in front of all dE counters. The distance between 
the dE-counter or E counter and the beam axis is 32.5 mm and 37.1 mm, respectively. These 
detectors are fixed to a cylindrical and four faced base of oxygen free high conductive copper 
(OFHC). It is helpful to avoid microphonic noise and natural background. 

F. Data acquisition system 

Analog signals from each detector are fed into preamplifiers (model 142IH, ORTEC, for dE 
counters model 142B, ORTEC, for E-counter) with inorganic coaxial cables (Cu-S5ESS-05 DIP- 
SOL CHEMICAL Co., Ltd.). Since a shorter distance between the detector and preamplifier is 
desirable to reduce electrical noise, the distance is 45 cm. The signals from the preamplifiers are 
amplified by both spectroscopy amplifiers (SAMR model 472 ORTEC) and timing filter ampli- 
fiers (TFA, model 454, ORTEC). The signals from the SAMP are led to a CAM AC peak sensitive 
analog-digital-converter system (ADC ,model AD811, ORTEC) while the signal from the TFA's 
are sent to a system of constant fraction discriminators (CFD, model 935, ORTEC) where the 
thresholds are set above the noise level of the detectors. The logic output of the CFD is fed into 
a logic Fan-In/Fan-Out (Fan-I/O, model 429A, ORTEC) delivering a gate signal for the CAMAC 



9 



ADCs with a gate width of 400 nsec. The logic output provides a start signal for a time spectrum 
via a CAMAC time-digital-converter system (TDC, model 2228A, Lecroy). The stop signal of 
the TDC's is provided by the CFD with a 100 nsec delay. The signals from the ADC and TDC 
systems are contlolled by a CAMAC crate controller (model CC7700, TOYO). The data from the 
crate controller is transferred to Linux station (model L400c, DELL), via a CAMAC bus, and 
stored on hard disk. A schematic diagram of the present data acquisition system is shown in Fig. 


The dead time of this data taking system is 400 /xsec for 1 event. The typical counting rate of 
the measurement of the 3 He + 3 He reaction, which includes the background events caused by 3 He 
+ d reaction, cosmic ray, and electrical noise is usually about 3 counts/sec. Therefore, the total 
dead time of these measurements is about 0.1 %. 

Before the reaction experiments, all the counters were calibrated using a 241 Am a source (5.48 
MeV). The energy resolution of the dE-counters was 100 ~ 120 keV (FWHM) for a 5.48 MeV 
a particle, and for the E-counter it was 70 ~ 100 keV(FWHM). The energy gain of S-amp was 
optimized to be able to measure the energy range up to 10 MeV for the dE-counter and 20 MeV 
for the E-counter. This energy range is required to measure not only 3 He + 3 He events but also 3 He 
+ d events, since 3 He + d events are needed for estimation of background and they are useful for 
checking the energy scale of each counter telescope. The linearity of the present amplifier system 
was measured by use of a precision pulse generator (model 419, ORTEC). Linearity was observed 
to within 0.05 % for all of the energy range of each counter. The stability of the energy gain of the 
amplifiers was checked and the resultant gain shifts were less than 2 % for 6 months. 

III. ANALYSIS OF EXPERIMENTS 

For evaluation of the cross section, the number of counts for the 3 He( 3 He,2p) 4 He reaction, the 
3 He target density and the 3 He beam intensity should be measured. The values for the effective 
reaction energy and the integral term for detection efficiency were calculated by means of a Monte 
Carlo simulation developed in the frame of the present work. 

The number of counts dN(z) per unit of time with respect to a length dz of the extended 3 He 
gas target is given by the expression 

dN{z) =N t -N h - a{E{z)) ■ rj(z) ■ dz, (3) 
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where N(z) is the number of counts for the 3 He( 3 He,2p) 4 He reaction, N t is the 3 He target density, 
iVb is the 3 He beam intensity per unit time, and 77(2) is the absolute detection efficiency. 

Introducing the stopping power e (i.e. the energy loss per unit length), equation © can be 
rewritten in the form 

dN(E) = N t -N h - a(E) ■ r](E) ■ e(E)- 1 ■ dE. (4) 

The total number of counts for the full target length is then given by 

N = N t -N h - j a(E) ■ V (E) ■ e(E)- 1 ■ dE. (5) 

For the case of a thin target, introducing an effective reaction energy E cS corresponding to the 
mean value of the projectile energy distribution in the detection setup, one arrives at 

N = N t -N h - a(E cS ) ■ J rj(E) ■ e(E)- 1 ■ dE. (6) 
A. Effective reaction energy 

The effective reaction energy was the mean value of the beam energy derived from an energy 
loss calculation in the target gas. As we could not apply any measurement for the absolute en- 
ergy of ion beam, such as time of flight techniques or an Wien filter, we determined the value by 
measuring a voltage divided with a precise register chain of a ratio of 1/10000(STANDARD EN- 
ERGY, S-100) for the applied voltage to the ion source (SPELLMAN, SL-1200 (60 kV/20mA)). 
This resistive voltage divider was investigated by applying the exact voltage calibrated with the 
second standard and the resultant absolute accuracy is ±0.2 %. It was measured at intervals of 1.5 
second for all measurements. The stability of the voltage was less than 0.1 % for about one day. 

Compared with an ion source of a quiescent plasma such as a duoplasmatron ion source an 
ECR ion source has a finite plasma potential. Thus, we took into account this plasma potential for 
the acceleration voltage. The adopted value was 21.3±2.4 eV as reported by JAERI group for the 
NANOGUN ECR ion source fl. 

For the extended geometry in the present gas target experiment, the reaction energy distribution 
due to the energy loss of the ion beam along the beam path should be estimated precisely as 
possible. In low energy experiments this might raise the ambiguity for the electron screening 
potential; we had to take care of experimental conditions such as target pressure or its difference 
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along the beam axis. In this series experiments for energies less than 30 keV center of mass energy, 
this problem should be treated more rigorously. 

There are quite number of experimental and theoretical works for stopping powers of charged 
particles in matter. Charged particles lose their energy through collisions with nuclei and with 
atomic electrons in matter. Although the greatest part of the energy loss occurs by collisions with 
electrons, low-energy ions lose their energy by collisions not only with electrons but with nuclei. 
Since it is impossible to deduce stopping powers data near the zero energy, using present-day 
technology, we have to use updated compilations with an accuracy of ranges between 2 % and 10 
Therefore to calculate the energy loss in the target, we used the stopping power values 
estimated by the SRIM computer code [11], which gives results which are consistent with the 
experimental energy losses to within a 10 % difference at most. 

For example, the stopping power of incident 3 He with an energy of Ei a b=90.00 keV is 
9.3 x 10~ 15 eV/atom/cm 2 . The energy distribution of this particle in the 3 He gas target with pres- 
sure of 0. 1 Torr is simulated by the Mote Carlo program and the result is shown in Fig. [10| 

We employ the full target length of L =30 cm as the distance between the entrance of the target 
just after the collimator and the entrance of the beam calorimeter. Since the rapid reduction of the 
cross section is about 1 1.2 % at E L AB( 3 He)=90 keV over the target thickness when we assumed a 
constant S(E) factor, the effective energy loss is evaluated to be (500±50) eV. This value roughly 
agrees with the back-of-the-envelope value of 493 eV for the energy loss between the entrance 
and the center of the counter telescope on the beam axis. As shown in Fig. |T0l we simulated 
the effective target length and estimated energy spread of 79 eV as the error of the incident beam 
energy. The energy loss due to the residual gas between the ion source and the target entrance is 
3.7 x 10~ 3 % of incident beam energy. 

In summary, at an incident beam energy of (90.0±0.13) keV, the effective reaction energy is 
Eiab = (89.50±0.13 ) keV, taking into account the accuracy of 0.1 % for the acceleration voltage, 
10 % for the stopping power and 0.09 % for the energy spread in the target. 



B. Beam current 



The incident projectile number (Nb) is deduced from the deposited power measured with our 
calorimetric device as described in the former section. It was calibrated using a charged beam in 
vacuum by comparing with the electrical beam current in the Faraday cup. The electrical charge 
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collected in the Faraday cup was measured with the current integrator (KEYTHLEY 616 digital 
electrometer). The absolute value of the current integrator was calibrated by measuring the current 
which was supplied with a precise current source (R6161 ADVANTEST). It has an accuracy and 
a stability better than 0.001 %. The difference between current measured by current integrator and 
the value of the current source is less than 3 % and it was corrected. The following was taken into 
account to evaluate the incident particle number Nf,. The beam current measured by calorimeter 
was corrected because of the energy loss of incident beam in the target gas. The energy loss was 
estimated by SRIM program and was 1.60±0.16 keV when incident beam energy was 90.0 keV 
(Ecm=45.0 keV). This energy loss was 1.8 % lower than the incident energy. 

The intensity was simultaneously corrected by recording the beam energy and the target gas 
pressure at intervals of 1.5 seconds. The typical beam intensity was about 100 fj,A at an incident 
energy of 90 keV for the 3 He 2+ beam. The beam intensity measured by the present system during 
experiment is shown in Fig. [TTJ 

C. Target density 

There are several factors which affect the target density N t such as the gas temperature and a 
pressure gradient in the target chamber. We measured the target gas temperature with a termister 
(103AT-2) inside the chamber which was likely to be different from that of the laboratory room 
since the target gas was heated by the beam and cooled by the circulated gas for the purification 
system. 

As the target pressure could not be measured directly at the beam-target interaction region 
during the experiment, the pressure was measured at the top of the target chamber as shown in 
Fig. [T21 The pressure distribution caused by the geometry of detector holder, collimators and gas 
circulation was measured by extending the stainless tube directly from the capacitance manometer 
set downstream of the target chamber. The capacitance manometer, which is usually installed at 
the top of the target chamber, was removed to the end of the chamber only at this measurement, as 
shown in Fig. [T2J Simultaneously, we used another gauge just before the inlet to the chamber for 
a normalization. The difference between the target gas pressure at the top of the chamber and that 
measured at the interaction region was rather small, which might be a shorter mean-free path at 
gas pressure of the order of 0.1 Torr. The absolute pressure at the target should be corrected by 5 
% less than that measured at the top of the chamber. Owing to these corrections, the target density 
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can be determined to an accuracy of 1.3 %, considering the accuracy of 0.16 % from the target gas 
pressure, 1 % from the correction due to the gas temperature, and 0.8 % from the correction due 
to the measurement position. The measured and corrected target gas pressure at the every interval 
of 1.5 seconds is shown in Fig. [T31 

IV. DATA ANALYSIS 

A. Monte Carlo simulation for the OCEAN experiment 

In order to find an optimum detector setup for a high efficiency and background free measure- 
ment, we exploit the Monte Carlo simulation program based on GEANT3. It was used to calculate 
the interaction between the ejectiles and the detectors. Also, the GENBOD code was used to 
generate the ejectiles. Thus the program takes into account the following aspects: 1) the detector 
geometry, 2) the energy loss and energy straggling of the ejectiles in both the target gas and i the 
thin foil in front of the detector, 3) kinematic effects on the energy of the ejectiles in the target, 4) 
yield dependence of the ejectiles over the passage of the target, and 5) the non-uniform depletion 
thickness for the E counter. 

B. Measurement of the D( 3 He,p) 4 He reaction 

To verify the validity of the simulation program, that is to estimate the systematic error, the 
experimental results of the D( 3 He, p) 4 He reaction (Q=18.4 MeV) are compared with those of 
simulations. The comparison to this reaction has several advantages: 1) The generated protons 
from the 3 He + d reaction have a definite energy of 14.7 MeV. 2) The energy of the protons from 
this reaction is almost the same as that from the 3 He + 3 He reaction. 3) the cross section of the 
3 He + d reaction is six orders of magnitude larger than that of 3 He + 3 He reaction. 

The D( 3 He,p) 4 He reaction was performed by using 90 keV (E cm =45 keV) 3 He2+ beam at the 
OCEAN facility. The target pressure of the deuterium gas was maintained around l.Ox 10~ 4 Torr. 

Fig. [J31 shows observed and simulated energy spectra obtained with the dE-E are compared. 
The broad energy spectra for E = 5~14 MeV at 5E = 1 MeV arise from an insufficient depletion 
depth for protons incident on the surface at angle near 90 degree. We applied a bias voltage of 180 
V to the E-counter, which corresponds to the depletion depth of 900 fim, to avoid a discharge in the 
gas target thicker depth needs more voltage. The other two structures arise from kinematics effects 
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in combination with protons which are incident at angular ranges of 135°~180° and 0°~45°. 
These three remarkable features are well simulated in the energy spectrum as shown in Fig. [14] 
(b). 

C. Background analysis 

It is crucial for the present measurement of the 3 He( 3 He,2p) 4 He reaction to identify the back- 
ground origin and to discriminate the true events from the fake events. As already stated, deuterium 
contamination in the target are the most serious. The number of deuterons in the gas target was 
determined from the data during the measurement of the 3 He( 3 He,2p) 4 He reaction as shown in 
Fig. [15] For this estimate, the value of the cross section for the D( 3 He,p) 4 He reaction was taken 
from ref. I12I1 . We conclude that the deuterium contamination is 0.2 ppm in the target gas, and 
that such a level could make background events at least 0.1 % of the observed events of the 3 He 
( 3 He,2p) 4 He reaction at the energy of E cm =45 keV. 

Another source of background events arises from electrical noise and cosmic rays. These are 
observed during the measurement without 3 He beam for 38 days of operation of OCEAN as shown 
in Fig. [16] The contribution from this background to the window of the 3 He( 3 He,2p)a; reaction 
is 3.6 counts/day. Of these, the cosmic muon events are located around dE = 70 keV and E =450 
keV because of the minimum ionization loss of 2 MeV cm 2 /g. We attempted to reject these events 
by applying the veto-counter upper and lower places for the target chamber. Finally the expected 
rate of the present reaction at lower energy is around a few events per day or less, and a typical 
single background rate of silicon detectors is one event per hour or more. In order to remove such 
accidental events, two proton coincidence should inevitably be required for the identification of 
the present reaction near the Gamow peak (next series of OCEAN experiment). 

D. Detector efficiency 

We developed a reasonable method to determine the acceptable area for the real events of 
the 3 He( 3 He,2p) 4 He reaction in the dE-E scatter plot without a redundant and ambiguous pro- 
cedure. Four types of data such as observed events of the reaction, simulated events for the 
3 He( 3 He,2p) 4 He reaction and for the D( 3 He,p) 4 He reaction, and observed background events are 
summarized for the analysis of each experimental run. The energy distribution of the dE-E scatter 
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plot is divided into 16000 parts of 100 keV x 100 keV parts as shown in Fig. [T7| The signal to 
noise ratio which should be derived from the Monte Carlo MC ( 3 He + 3 He ) divided by MC ( 3 He 
+ D) and measured background was allotted for 16000 parts. All parts are ordered as a function of 
their S/N; parts having a better S/N are located at the right hand side while worse parts are located 
to the left, as shown in Fig. [TSJ Fig fTST a) shows the distribution of the simulated events from the 
3 He + 3 He reaction as a function of S/N ratio. Many events are located in the right hand side in the 
figure, that should correspond to the better S/N. Also, Figs. \T%\ b) and (c) show the distribution of 
the simulated events for 3 He + d reaction events and the observed background events as discussed 
above. It is usual that these background events should be located in the left hand side of each 
figure. Fig. [TST d) shows the distribution of the observed event for 3 He + 3 He reaction at E cm =45 
keV, as a function of S/N. The contribution from the background events is apparently very small. 
Therefore the observed distribution as shown in Fig. fTSl' d') is particularly similar to the simulated 
one, shownn in Fig. HHta), without subtraction of any background events as shown in Fig. [T8l cD. 

Since, most of the background events exist in channels less than 13000 (Fig. IT8l the acceptable 
area for the 3 He + 3 He reaction could be assigned to channels larger than 13000 channel. The 
region of that is shown in Fig. fT9l 

Experimental results are as follows; 3344 counts are observed in the acceptable region, while 
the contribution from the 3 He + D events to the region is 20.9 counts, and that from the other 
background component is 2.46 counts. After subtracting the number of these background events 
from the number of observed events in the acceptable region, the number of true events for the 
3 He + 3 He reaction is 3337.4, with a statistical error of 1.8 %. According to this procedure, the 
detection efficiency r](x) can be written as 

16000 

V (x) = (7) 

i=x 

where i = x to 16000, N & {i) is the number of counts for the simulated distribution, and x is the 
parameter of the boundary cut point for the accepted events. The accepted event have been derived 
as 

16000 

A(x) = - N h (i) - N c (i)}, (8) 

i=x 

where, A(x) refers to accepted events for the 3 He( 3 He,2p) 4 He reaction, Nd(i), Nb(i), N c (i) are 
the number of counts for the observed events, for the d( 3 He,p) 4 He and for other background 
events, respectively. Finally the ratio of A(x)/y(x) corresponding to the cross section of the 
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3 He( 3 He,2p) 4 He reaction can be obtained. The ratio slightly depends on the boundary parameter 
x as shown in Fig. |20| Thus, the accuracy of the simulated energy distribution in a scatter plot 
should be derived from the fluctuation of this ratio. When we include the geometrical uncertainty 
of the counter telescope, i.e., 0.5 % uncertainty of the detection efficiency simulated with Monte 
Carlo program, the systematical error of the detection efficiency is evaluated to be 3 %. 



V. EXPERIMENTAL RESULTS AND DISCUSSION 



The experimental results for the cross section and S-factors obtained for the first series of 
OCEAN experiments as from the year 2000, together with the experimental conditions such as 
the live time, beam current, the target gas pressure, and the target temperature are shown in table 
Hill The observed events for the 3 He( 3 He,2p) 4 He reaction and the background events from various 
sources are also shown in the same table UTTl 

The cross section for the 3 He( 3 He,2p) 4 He reaction has been derived from the following equa- 
tion, 

N = N t -N h - a(E cS ) ■ J r](E) ■ e(E)- 1 ■ dE. (9) 

where N is the number of counts for 3 He( 3 He,2p) 4 He reaction, N t is the 3 He target density, N h 
is the 3 He beam intensity, E c s is the effective reaction energy, i](E) is the absolute detection 
efficiency, and e(E) is the stopping power. 

The astrophysical S-factors were deduced from the equation 

a(E) = ^ex P (-2n V ), (10) 

where r] is the Sommerfeld parameter given by: 

27T77 = 31.29Z 1 Z 2 (^) 1 / 2 (11) 
E 

where Z\ and Z 2 are the nuclear charges of the interacting particles in the entrance channel, [i is 
reduced mass (in units of amu) and E is the center of mass energy (in units of keV). 

Our present data between E cm = 45 keV to 3 1 keV of S(E) are in good agreement with the results 
of existing data of Krauss et al The accuracy of both statistical and the systematical data of 
the present measurement is better than that in ref. |Q]. 

For the last two decades, studies of the 3 He( 3 He,2p) 4 He reaction have been carried out over a 
wide range of energies (Bacher and Tombrello 111 311 - Dwarakanath and Winkler Krauss et al 
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and LUNA [1]). Of the previous studies, experiments between 17.9 to 342.5 keV center of 
mass energy by Krauss et al cover a wider energy range than others since they used two accel- 
erators of 350 kV accelerator in Munster University and 100 kV facility at Bochum University. 
Recently the LUNA group in the LNGS has presented data down to 16.50 keV center of mass 
energy, although the data was obtained with two separate experimental set-ups,namely at the 450 
kV accelerator at Bochum and also 50 kV at LNGS. There remains unaccesible region of energy 
between 45.82 keV and 24.80 keV. The data in this energy region has been supplied by Krauss et 
al. 

In these previous experiments there were many efforts to obtain the true scarce events from the 
contaminated background events, such as deuterium contamination both in the incident beam and 
in the target gas, cosmic rays (mainly muon) or heavy particles, and electric noise. There have 
been several solutions for these difficulties; Krauss et al., pointed out in 1987 the purity of the ion 
beam and of the target were of special interest and they estimated that at 350 keV the mass-3 beam 
contamination HD + was of the order of 10 5 [2]. Also they applied proton-proton coincidences to 
discriminate the real events from intruded events for the D( 3 He,p) 4 He reaction. They surrounded 
the target chamber with NE102A plastic scintillator in order to identify for cosmic events. In this 
way the unidentified cosmic background coincidence event rate was estimated to be less than 1 
event/200h in the measurements at E cm < 50 keV. Another source of background events due to 
occasional discharges of the accelerator high voltage could be avoided by amplifier filter or noise 
filter. In spite of these efforts, they still consider the contribution of the events caused by the high 
energy protons from the 3 He(d,p) 4 He reaction (e.g. 2.5 % contribution at E cm =25 keV ). They 
also evaluated that the background contribution to the 3 He+ 3 He region is 0.40 % of the observed 
counts of the d- 3 He events. 

The measurements at the DTL showed that cosmic rays produced events within the dE-E region 
of the 3 He( 3 He,2p) 4 He reaction amounts to 3.5 x 10~ 4 events/s, while in the latter place at LNGS 
this rate was observed to be reduced by at least a factor 200, that is, 1 .8 x 10~ 6 . It will be negligibly 
small for the low energy measurement less than 30 keV center of mass energy. 

We now argue that the OCEAN facility could overcome these difficulties by applying the fol- 
lowing, while there still exist the possibilities of background due to the target gas: 

1) Only OCEAN exploits doubly charged He-3 ions produced with an ECR ion source and it 
could avoid the background events due to D (beam) + 3 He. Krauss et al performed a p-p coinci- 
dence experiment to discriminate the true event, but it could not be applied for the center of mass 
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energy less than 25 keV. 

2) Our facility OCEAN has been located in an experimental area of the cyclotron and shielded 
with 5 m of thick concrete. We have observed the background events arising from cosmic rays in 
38 days and also observed the effect from the induced events due to the high energy particles by 
using Veto-counters composed of plastic scintillators located up and down the target chamber. By 
this condition of the OCEAN sight we still have a sufficient signal to noise ratio of 10 times larger 
than unknown fake events as listed in table Hill 

Generally, there are still many discussions for the screening potential which enhance the cross 
section remarkably for low energy fusion reactions. A recent result of the 3 He( 3 He,2p) 4 He reaction 
data assumed for electron screening with U e = 330 eV and adopted a polynomial function for S(E) 
factor,namely S(E) = 5.18 + 2.22E + 0.804E 2 . The S(0) value is close to the recommended value 
of AD98 (5.4 + 0.4 MeV b). In a review article concerning synthesis of the elements in stars 
it cited LUNA data as S(E ) = 5.3±0.05(stat.) ±0.30(syst.) +0.30(C/ e ) MeVb. It still includes 
a U e value which suggests that the effective value of U e may be larger than the adiabatic limit. 

Furthermore, it is mentioned that improved data should be obtained at energies from E = 25 
keV to 60 keV to confirm or reject the possibility of a relatively larger systematic error in the S(E) 
data around these energies. 

VI. CONCLUSION AND PERSPECTIVES 

The present experiment with OCEAN has proved to reinforce the compilation data for nuclear 
astrophysics in solar fusion rates. High current and low background system could measure the 
cross section in the center of mass energy range from 50 to 30 keV. As a second series of experi- 
ments, we have started to measure at energies less than 30 keV, and these results will be reported 
in the near future. Better results with respect to reduced stematic and statistical errors compared 
to existing data are expected. When we apply a beam intensity of 1 mA for the 3 He 2+ beam at E cm 
= 20 to 30 keV and a target gas pressure of 0.1 Torr, we can expect 70 ~ 2 events per day for real 
events, while the fake event will be 6 ~ 1.8 events per day. Hence, we expect a measurement with 
a statistical error of about 10 % in this energy range. 
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TABLE I: The achieved beam current at the ion source and the target. 



H.V. 


Target(^A) Ion source(/iA) 


40 kV 1+ 


1208 


3010 


2+ 


103 


3000 


30 kV 1+ 


1200 


3800 


2+ 


35 


1900 
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TABLE II: Comparison of the heat calculated from the kinetic energy of incident particles and from the 
HFS output. The beam current is measured by means of a Faraday cup, and the Voltage at the IS(Ion 
source) determines the kinetic energy of incident particles. The beam power is calculated by multiplying 
these factors. The heat through the HFS is calculated using the calibration parameter, which is given as 
1.10~1.11 /iV/W/m 2 @70°F 



Beam Current Voltage at IS HFS Output Energy of Beam Heat Transfer 



(//A) 


(kV) 


(mV) 


(W) 


(W) 


0. 


0. 


-0.83 


0. 


0.17 


96.4 


40.0 


20.86 


3.84 


3.75 


69.7 


20.0 


7.65 


1.39 


1.37 


74.3 


25.0 


9.88 


1.86 


1.78 


82.5 


35.0 


15.66 


2.89 


2.81 
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TABLE III: Summary of measurements for the 3 He( 3 He,2p) 4 He reaction. L.T. : Live Time, B.C. : Beam 
Current, T.P. : Target Pressure, T.T. : Target Temperature, C.S. : Cross Section, True : 3 He+ 3 He, BGl : BG 
3 He+d, BG2 : BG Other, Cnt : Counts, S-Fac: Astrophysical S-factor. 





L.T. 


B.C. 


T.P. 




T.T. 


True 


BGl 


BG2 


C.S. 




S-Fac 


(keV) 


(sec) 


QjlA) 


(Torr) 




(°C) 


(Cnt) 


(Cnt) 


(Cnt) 


(barn) 




(MeV-b) 


45.3 


92567 


104. 


7.46x10" 


-2 


27.1 


3276 


20.9 


2.46 


1.53x10- 


-8 


5.39±0.09 


43.3 


78647 


91.4 


6.72x10- 


-2 


27.3 


1374 


7.50 


2.09 


9.55x10" 


-9 


5.43±0.14 


41.3 


80687 


100. 


6.74x10- 


-2 


27.1 


939 


7.08 


2.15 


5.79x10- 


-9 


5.51±0.18 


39.3 


83109 


87.4 


7.26x10" 


-2 


27.0 


542 


6.08 


2.21 


3.44x10- 


-9 


5.69±0.25 


37.3 


155442 


112. 


8.24x10- 


-2 


29.3 


770 


17.0 


4.14 


1.83x10- 


-9 


5.46±0.20 


35.2 


338862 


100. 


8.21x10- 


-2 


29.3 


770 


21.4 


9.02 


9.46x10- 


-10 


5.62±0.21 


33.1 


615814 


103. 


8.25x10- 


-2 


30.4 


691 


11.4 


16.4 


4.52x10- 


-10 


5.48±0.22 


31.2 


528134 


93.6 


8.28x10- 


-2 


30.3 


293 


5.02 


14.1 


2.46x10- 


-10 


6.40±0.39 
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FIG. 1: Complete layout of Osaka University Cosmological Experimental Apparatus for Nuclear Physics, 
OCEAN. 
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FIG. 2: Beam transport calculation for the present designed scheme consisting of a 90 ° dipole mag- 
net+quadrupole doublet for 1 mA 50 keV 3 He 1+ beam. The beam size at the target (vertical and horizontal) 
was calculated as a function of the distance between the ion source to the dipole magnet or dipole magnet 
to Q-magnet, respectively. The beam source is assumed to be 100 7rmm mrad and 5 mm in diameter. 
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FIG. 3: Complete layout of the window-less gas target evacuation and recirculation system. TMP: Turbo 
molecular pump. R.P: Rotary pump. 
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FIG. 4: Cross sectional view of the calorimeter. 
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FIG. 5: (a)Measured beam current(/(^A)) v.s. heat Hu\(H (mV)), under the condition of SH/St = in 
eq. (2). The parameter k\ in eq. (2) is derived by fitting with a linear function. (b)Plot of SH/5t(eq. (2) in 
text) as the difference of / — k\ ■ H. The parameter /c2 in eq. (2) is derived by fitting with a linear function. 
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FIG. 6: Beam current, Q) measured by the Faraday cup, Q, calculated from k\-H and Q>ki-H+k2-5H / 5t as 
a function of time. The measurements were carried out for a)stable beam current b)increasing beam current 
c)decreasing beam current and d)for rapid beam step of current. 
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FIG. 7: The accuracy of the beam current for different beam energies 35, 30, 25, and 20 keV, in the form 
of (IHFS-IFC)/IFC, where IHFS and IFC denote the beam currents measured by the HFS and the Faraday 
cup(FC). 
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FIG. 8: Schematic view of detector assembly. 
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FIG. 9: View of electronics for data acquisition. 
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FIG. 10: Simulated interaction energy distribution for 3 He+ 3 He reaction at Ei a t>=90.0 keV. A target gas 
pressure of 0. 1 Torr is assumed. 
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FIG. 11: Measured and corrected intensity of 3 He 2+ beam at Ej a b=90.0 keV. Measurements were made at 
1.5 second intervals. 



35 



Target gas 



Detector 



Detector support 




Target chamber 
Interaction region 



Capacitance 
manometer 



i i i i i i 

10 cm 



FIG. 12: Schematic view of the measurement of the gas pressure at the interaction region between the beam 
and gas target. 
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FIG. 13: Measured target 3 He gas pressure by Baratron capacitance manometer. Measurement intervals 
were 1 .5 sec for all experiments for the 3 He+ 3 He reaction. The pressure was normalized to that of temper- 
ature at °C. 
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FIG. 14: The observed (left) and simulated (right) energy spectrum for the d( 3 He,p) 4 He reaction at 
E3# e =90.0 keV by means of the AE-E counter telescope. 
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FIG. 15: Observed AE-E energy spectrum of the 3 He( 3 He,2p) 4 He and D( 3 He,p) 4 He reactions. The solid 
line shows the acceptable energy region for the D( 3 He,p) 4 He reaction. The deuteron contamination in the 
target 3 He gas was evaluated from the events in this region. 
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FIG. 16: Background events arising from electronic noise and cosmic-ray, observed with the same geometry 
as the 3 He( 3 He,2p) 4 He experiment. 



40 




8 10 12 14 16 
E-counter energy(MeV) 



FIG. 17: Schematic view of the the acceptance region. The AE-E energy scatter region was divided into 
16000 parts of 100 keV x 100 keV divisions. Signal to noise ratio was examined in each part. 
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FIG. 18: Event distribution ordered as a function of their S/N ratios for four types of data. The parts having 
a better S/N locate at the right hand side while worse parts are located to the left, (a)simulated 3 He+ 3 He, 
(b)simulated 3 He+d, (c)other observed background(electric noise and cosmic-ray), (d)observed 3 He+ 3 He 
are shown. S/N is given by Stimulated 3 He+ 3 He)/ N(simulated 3 He+d + observed other background). 
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FIG. 19: AE-E scatter plot obtained from the 3 He+ 3 He reaction experiment. The region assembly of 
grosspoints shows the accepted region as the true events from the 3 He( 3 He,2p) 4 He reaction. 
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FIG. 20: Ratio of A(x)/r)(x) as the function of the boundary parameter x. The ratio corresponds to 
the cross section of the 3 He+ 3 He reaction. r](x) and A(x) was derived from eqs. Q and (JSJl [see text], 
respectively. 
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